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Abstract 

We propose a new mean-field-type framework which can treat the strong correla- 
tion induced by the tensor force. To treat the tensor correlation we break the charge 
and parity symmetries of a single-particle state and restore these symmetries of the 
total system by the projection method. We perform the charge and parity pro- 
jections before variation and obtain a Hartree-Fock-like equation, which is solved 
self-consistently. We apply the Hartree-Fock-like equation to the alpha particle and 
find that by breaking the parity and charge symmetries, the correlation induced by 
the tensor force is obtained in the projected mean- field framework. We emphasize 
that the projection before the variation is important to pick up the tensor correlation 
in the present framework. 
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1 Introduction 

Mean-field models witli effective interactions like the Skyrme force [1] or the 
Gogny force [2], are widely applied to various fields in nuclear physics with 
great success. The relativistic mean field theory is also applied to many nuclear 
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systems successfully [3,4]. These models are based on the single-particle picture 
of the nuclear many-body system. In this picture, a nucleon in a nucleus is 
moving in a nuclear mean field made by other nucleons in the nucleus. The 
single-particle picture is known to be appropriate in the wide range of nuclear 
phenomena since the estabhshment of the shell model [5] . 

Usually, the effective interactions used in nuclear mean-field models only have 
the central and LS parts as the nuclear two-body interactions. They do not 
have the tensor part. However, the tensor force is known to be very important 
in the structure of nuclei [6,7,8,9,10,11,12]. The tensor force is mainly mediated 
by the pion, which is one of the most important meson since Yukawa proposed 
it as a mediator of the nuclear force [13]. Many studies have demonstrated that 
the tensor force gives large attractive energy in light nuclei [6,7,8,9,10]. The 
shell model calculation shows that the tensor force can produce about a half of 
single-particle Is splittings in light nuclei [11]. These facts show that the tensor 
force is an essential ingredient of the structure of light nuclei. Furthermore, the 
tensor force plays an important role in the saturation mechanism of nuclear 
matter [12]. This fact indicates that the tensor force may have an important 
effect also on the structure of nuclei in the heavy mass region. 

The importance of the tensor force mentioned above inspires us to treat the 
tensor force seriously even in a mean-field model. Many people, however, be- 
lieve that the effect of the tensor force can be treated by renormalizing the 
central and LS parts of the effective two-body interaction. This renormaliza- 
tion is founded by the G-matrix theory, where the effective interaction is made 
by renormalizing the central and LS forces to include the effect of the tensor 
force together with the short-range correlation. In the G-matrix theory, the 
effective tensor force is weakened and treated simply as the residual interac- 
tion. The main part of the tensor force is included in the '^E part of the central 
force. We see the importance of the tensor force by comparing the theoretical 
results for nuclear matter and the alpha particle. In nuclear matter, the ratio 
of the energy contribution from the part of the effective interaction and 
that from ^E is nearly equal to 1, but in the alpha particle, this ratio increases 
to about 1.5 [6]. This enhancement of the energy contribution from the ^E 
component indicates that the effect of the tensor force is larger in the alpha 
particle than in nuclear matter, and the tensor force should be more important 
in finite nuclei than in nuclear matter. 

There are some Hartree-Fock calculations using the nuclear forces which in- 
clude the tensor force [14,15,16] or the pion [17,18] exphcitly. They show that 
the tensor force gives a large contribution to the Is splitting of single-particle 
states especially in the spin-unsaturated nuclei. In contrast to the conventional 
mean-field calculations, we recently proposed a new method to treat the pion 
in the relativistic mean field (RMF) theory [19]. Because the pion is a pseu- 
doscalar meson, a single-particle state of a nucleon in a nuclear mean-field 
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should change its parity when the nucleon emits or absorbs a pion. It suggests 
that we can incorporate the pion into the RMF theory by introducing parity- 
mixed single-particle states. The idea of the parity-mixed single-particle state 
was already discussed in the references [20,21,22,23,24] many years ago. In the 
previous paper [19] , we apphed this idea in the RMF theory and found that 
the pion mean field becomes finite, due to which the correlation induced by 
the pion is appropriately taken into account. This result is not restricted to 
a relativistic framework. We can treat the tensor force by the introduction of 
the parity-mixed single-particle states in a non- relativistic framework also. 

The purpose of this paper is to apply the same idea in a non-relativistic mean- 
field model and see the effect of the tensor force in the alpha particle. We make 
two major improvements compared to the previous study in the treatment of a 
nuclear many-body system. The first improvement is straightforward as com- 
ing from the isovector character of the pion. Because the pion is a isovector 
meson, a single-particle state in a mean field changes its charge state when it 
emits or absorbs a charged pion. This fact indicates that the isovector nature 
of the pion can be treated by introducing charge-mixed single-particle states, 
which have both proton and neutron components. The dominant part of the 
tensor force is Ti • T2-type and therefore this improvement is important when 
we treat the tensor force. The second improvement is related with the restora- 
tion of symmetry. We mix the parities and the charges into a single-particle 
state. A total wave function made from the single-particle states with the 
parity and charge mixings does not have a good parity and a definite charge 
number. Therefore we need to perform the parity projection and the charge 
(number) projection to obtain a total wave function having a good parity and 
a definite charge number. By performing the projections we can pick up the 
2-particlc-2-hole, the 4-particle-4-hole, . . . correlations, which are induced by 
the tensor force for the ground state of even-even mass nuclei [19], and are 
important in the binding mechanism from light nuclei to nuclear matter. 

We may think of performing the charge and parity projections after we obtain 
a total wave function assuming the charge and parity mixings (projection 
after variation). Because the effect of the tensor force is strong and affect 
the binding mechanism of nuclei largely, we need to treat the restoration of 
the parity and charge symmetries carefully. Therefore, in a more sophisticated 
way, we should first perform the parity and charge projections, and then take a 
variation of the energy expectation value which is evaluated with the projected 
total wave function (variation after projection (VAP)). In the present study 
we perform calculations for the alpha particle in the VAP scheme and show 
that the VAP scheme is needed to treat the tensor force more appropriately. 
We should note that the present framework is not a simple mean-field model 
because we effectively superpose many Slater determinants by performing the 
projections. In this sense, our framework is a kind of descriptions beyond the 
mean field [25,26,27,28,29]. 
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This paper is arranged as follows. In Section 2, we will formulate the charge- 
and parity-projected Hartree-Fock method. In Section 3 we will apply the 
charge- and parity-projected Hartree-Fock method to the alpha particle to 
show the effectiveness of our method. In the last section we will give the 
summary. 



2 Charge- and peirity-projected Hcirtree-Fock method 

We present the charge- and parity-projected Hartree-Fock (CPPHF) method. 
In the CPPHF method, we construct a total wave function with single-particle 
states with the parity and charge mixings. A single-particle wave function with 
the charge and parity mixings in the spherical case is written as, 

t,=±l/2 

This wave function consists of four terms, each of which has a positive or 
negative parity, and a charge number, 1 or (proton or neutron). Here, the 
radial wave functions 0's depend only on the radial coordinate r, the isospin 
wave functions are denoted as C{tz) {tz = 1/2 for proton and tz = —1/2 for 
neutron), and yjim is the eigenfunction of total spin j = l + s. In the spherical 
case with the parity and charge mixings, the good quantum numbers are j 
and m. Two angular wave functions yjim and yfim are related to each other 
as 

yjiM = <T-ry,iU^), / = ||^| (i = !-J)- 

Here, cr is the Pauli spin operator and r is the unit radial vector. Because 
(T • f is a 0~ operator, yjim and yjim have the same total angular momentum 
j but different orbital angular momenta, i.e., \l — l\ ~ 1- It means that they 
have the opposite parities to each other. The symbol n is introduced to label 
the single-particle states with the same j and m. In the CPPHF method, n 
does not correspond to the node quantum number because of the parity and 
charge mixings. We take a Slater determinant made from the charge- and 
parity-mixed single-particle states as an intrinsic wave function for a nucleus 
with the mass number A: 

*^"*^ = ^^n^a„(x„). (3) 
V^! a=l 
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Here, ctfl denotes n, j, and m. in (1) and A is the antisymmetrization operator. 
As aheady mentioned above, ^I'™*'' does not have a good parity and a definite 
charge number. We perform the projections of parity (±) and charge number 
(Z) on 

^{±;Z) ^ pP(±)p'=(Z)*'"^*^ (4) 

Here, 'PP(±) is the parity-projection operator, where V^{+) projects out the 

positive parity state and V^{—) projects out the negative parity one. V'^{Z) is 
the charge-number-projection operator, which projects out the wave function 
with a charge number Z. Therefore, has a good parity (±) and a definite 

charge number {Z). The parity projection operator ^^(i) is defined as 

l±P ( ^ ^ \ 

rn±) = ^ [p=npaj, (5) 

where the total parity operator P is the product of the parity operator pa for 
each single-particle state. The charge projection operator P^{Z) is defined as 

V\Z) = — / dQe'^^-^^' = — / dQe-'^'C{Q) [z = y , (6) 

27r JQ 27r \ ^ 2 ) 

where Z is the charge number operator, which is the sum of the single-particle 
proton projection operator (l+rf ) /2, and the charge- phase operator is defined 
as C{B) = e'^^. 

We take a Hamiltonian H in the following form. 



A 

+ {^ck^ah) + ^^T(2;a6) + i)Ls(2;a6) + Vc^^x\{^ah)) (7) 

a>b=l 
j4 a 

= E^^W+ E ^M- (8) 

a=l a>b=l 

Here, the first and the second terms on the right-hand side are the single- 
particle kinetic energy and the energy of the center of mass motion. The 
two-body interactions, -Oc, "Ot, "^ls; and -Ocoui represent the central, tensor, LS, 
and Coulomb interactions, respectively. We define the one-body kinetic energy 
operator t{xa) and the two-body potential energy operator v{xah) for the later 
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convenience, 

, A-1 ^ 

v{Xab) = Vc{Xab) + VT{Xab) + VLs{Xab) + VConl{Xab) + o. " Vfo. (10) 

The energy correction due to the center of mass motion in (7) is included in 
i{xa) and v{xab)- 

We take the expectation value for a Hamiltonian H with the projected wave 
function and obtain the energy functional, 



/intr\ 



47r 



^ /o'"(i^e-^^^(n(o)(^)±n(P)(^)) 



(11) 



The denominator in the right-hand side of the above equation is the normal- 
ization of the total wave function, 

^(±;z) = ^^(±;z)|^(±;z)^ ^ ±_ /'"rf^e-^e fn(°)(^) ±n(P)(^)) . (12) 

Here, n'^'^^O) is the determinant of the norm matrix between the original wave 
functions ^-nd the charge-rotated wave functions ipa^(9). n^^\9) is the 
determinant of the norm matrix between the original wave functions ipaa 
the parity-inverted and charge- rotated wave functions ipaai^)- 

(13) 

The charge- rotated wave function ip^a i^b', d) and the parity- inverted and charge- 
rotated wave function i/j^^(xb; 9) are defined as 

,,^^{x,-ff)^^^^^)l^i,^^{xb), (14) 
V;£)(a;,;^)^p,e^^(^+-^^)/Va„(a;.), (15) 

where is the single-particle parity operator in (5) and (1 -|- t^)/2 is the 
single-particle proton projection operator in (6). 

The numerator in the right-hand side of (11) is the unnormalized total energy, 
(^(±;^)|^|^(±;^)) = i_^'"d^e-'^^ {E^^\e) ± E^''\e)) . (16) 
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E^^\9) in the right-hand side of (16) has a similar form as a simple Hartree- 
Fock energy but the single-particle wave functions in the ket are modified by 
the charge rotation, 

A A 
a=l a>b=l 

Here, ipaai^l ^) is the superposition of ipaai^] ^) weighted by the inverse of the 
charge-rotated norm matrix (i?(°)(6')~^)ba, 

^aSx\0) = J2^a,{^;0){B^'\e)-'\a. (18) 

6=1 

This summation for ip^^ [x; 9) comes from the antisymmetrization of the total 
wave function. E^^\d = 0) reduces to a simple Hartree-Fock energy, 
in the right-hand side of (16) has a similar form as E^^\9) but ^pa^{9ys are 
replaced by ip^^{9ys^ 

E^^\9) = \H PC {9)\¥'"'') 

-i:{i^aMSm+ E (V'a>aJi)|^S(^)^i?(^)-^SH^)^£n^)). (19) 
a=l a>6=l 

Here, ip^^{x] 9) is the sum of tjj^^^{x; 9) weighted by the inverse of the parity- 
inverted and charge-rotated norm matrix {B^^\9)~^)i,a, 

= j:^^^i{^;0){B^^\e)-'),a. (20) 

6=1 



We then take the variation of E^^'^^ with respect to a single-particle wave 
function ipaa, 

jj-- \e^^''^ - E ^^c{^Pa,\^aA = 0. (21) 
Oi/jkaiXa) [ b,c=l J 

The Lagrange multiplier tab is introduce to guarantee the ortho-normalization 
of a single-particle wave function, {il'aa\'4'at,) = ^aa,ab- As the result, we obtain 
the following Hartree-Fock-like equation with the charge and parity projec- 
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tions (the CPPHF equation) for each ipa^ , 
1 r^'^ 



Att Jo 



dOe 



-ize 



n 



(0) 



6=1 



^{ll^a,\v{Xal)\'l/jaa{^))l'^ab{Xa;0) 
6=1 

A 

Ml 



6=1 J 

I 6=1 

6=1 

6=1 

= n^^-^^^Y.^ab^-M, (22) 

6=1 

where a = 1,2,. . .,A. Here, rj^^i^) defined as follows, 



A 



+ T.{^aMv\^aM^aXO)-^aM^aM), (23) 



c=l 



+ E(^aaV'a.|i)|^£H^)v;i^ (^) - i^SWSi^))^ (24) 

c=l 

and the notation for the integration of the two-body matrix elements, 

{'ipab\HXal)\lpac)l = j dXi1pl^{xi)v{Xal)lpao{xi) ■ (25) 

The system of the coupled equations (22) for a = 1, • • ■ , A is solved self- 
consistently. We note here that the CPPHF equation reduces to the parity- 
projected Hartree-Fock equation with only the parity projection by setting 
^ = in (22), which was already obtained by S. Takami et al. [25]. 

We give here the expressions for the expectation value of the kinetic energy 
^^]^ ^jjg center of mass correction and that of the two-body potential 
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energy (va-)^^'^^ for ■Og. (cr = C, T, LS, and Coul) for the later convenience. 

^°(^)|i:(V'aJt|^a„(^)) 



1 r27c 



. a=l 



fc2 



a>6=l 



a=l 



a>b=l 



(26) 



X < n 



a>6=l 

±n(P)(^) E (^a>aJi).|^£H^)^i^W-^i^W^^^^ (27) 

a>&=l 



3 Application to the alpha particle 



We apply the CPPHF method to the alpha particle (A=4, Z=2), which has 
been studied extensively with various theoretical methods including exact cal- 
culations. As the ground state configuration of the alpha particle, we assume 
four states n = 1, 2 and (jm) = (1/2, ibl/2), which are fully occupied. Here, j 
and m are the total angular momentum and the magnetic quantum number 
of a single-particle state. The index n labels the two different states which 
have the same j and m. The intrinsic wave function (3) for the alpha particle 
becomes 

n n K=im{^)- (28) 

n=l,2 m=±i 

Because we fill m=±l/2 states with j = 1/2 for each n, "i/^"^^ has the total 
angular momentum 0. For each n, ipnj=^m becomes 



1 

2 

+<l>ni=uAr)yj=ii=immiQ). (29) 
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The single-particle wave function t/'^. i^^ consists of four components, which 
have different parities (±) and different charges (t^ = ±1/2). The first term 
on the right-hand-side has j = 1/2 and / = (a s-state with positive-parity) 
and the second term has j = 1/2 and / = 1 (a p-state with negative-parity). 
The single-particle p-state probability -P(— ), which is a measure for how much 
the p-state mixes into the simple (Os)^ configuration in the intrinsic state of 
the alpha particle, is defined as 




The factor 1/2 in front of the right hand side of (30) is added to average 
the single-particle p-state probability over the states for n = 1,2. The four- 
body intrinsic wave function (28) does not have a good parity and a definite 
charge number; it is a mixture of positive and negative parities, and of various 
charge numbers Z=0~4. In the alpha particle, the ground state has a positive 
parity and the charge number Z=2. Therefore, we act the projection operator 
P^{+)P^{2) on the intrinsic wave function (28). 

In the numerical calculation, we expand the radial wave function in (29) by 
the Gaussian basis with the geometric-series widths [30]. We take the number 
of basis as 10 and set the minimum width to 0.5 fm and the maximum width 
to 6.0 fm. To solve the CPPHF equation (22) self-consistently we use the 
gradient method [25,31]. 

For the central force, we take the Volkov force No. 1 [32] as a reference and 
introduce the multiplying factor xte for the triplet-even part. 

rj-{xa,) = - ^Aexp(-(r„,/aA)') (xteP^? + + (1 - 2mv)(P„? + P^)) 

+ vneM-(ra,/anr) {P'l + + (1 - 2mv)(P„? + P^^)) . 

(31) 

Here, Paf^ ^iFy P'ab'y Pah' ^^c projcctiou operators for the triplet- 
even, singlet-even, triplet-odd, and singlet-odd states, respectively. If we set 
xte = 1, ■^c™ reduces to the original Volkov force. The term proportional to vx 
is a middle-range attractive part and that to is a short-range repulsive one. 
The values of the parameters are = 83.34 MeV, «a = 1-6 fm, = 144.86 
MeV, and = 0.82 fm. We fix the Majorana parameter my = 0.6. The 
Volkov force reproduces the binding energy of the alpha particle without non- 
central forces. We can consider that the effect of the tensor force is included 
in the Volkov force by renormalizing the central force. This renormalization 
affects the middle range of the part of the central force mainly. Therefore, 
we reduce the "^E part of the central force in the Volkov force by multiplying 
the factor xte to the attractive ^E part of the Volkov interaction, when we 
include the tensor force explicitly. 
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As for the non-central forces (the tensor and the LS forces), the GSRS [33] force 
is adopted as a reference, which is determined so as to reproduce the scattering 
phase shift of NN scattering and have the one-pion-exchange potential tail in 
the tensor part. In the mean- field-type framework, we have to treat the effect 
of the short-range correlation on the tensor interaction and further take into 
account the effect of the A isobar- hole excitations [34], separately from the 
mean-field correlation. In the present study, we shall treat these effects by 
multiplying a factor to the Ti • T2 part of the tensor force, which is the 
dominant part of the tensor force and mediated mainly by the pion. We need 
further studies on this factor in the near future. The tensor part of the GSRS 
force with the multiplying factor can be written as following. 



3 3 



E vZ eM-irab/VTir) + 3 E < eM-irabhTn)') 

n=l n=l 
n=l n=l 



'-^abi 

(32) 



where the tensor operator STf^ is defined as. 



= 3 2 CTa-CTb- \o6) 



ah 



The factors, xte and xt are correlated and are determined so as to reproduce 
the binding energy of the alpha particle. 

The LS part of the G3RS force is made from the two-range Gaussians: 

^Ls(x„,) = E 4?„ ^M-{ro,hl%J)P'''La, . Sab. (34) 



n=l 



The parameters of the G3RS force used in this work is summarized in Table 1. 
We take further the Goulomb interaction, 

VCon\{x ab) = , (35) 

rab 

with e being the charge of the electron. 

We present here the results of the GPPHF calculations for the alpha particle. 
To see the effect of the charge and parity mixings, we calculate three cases. The 
first calculation is the simple Hartree-Fock (HP) scheme, in which we do not 
perform neither the parity nor charge projection. The second calculation is the 
parity-projected Hartree-Fock (PPHF) scheme, in which we only perform the 
parity projection. The third calculation is the charge- and parity-projected 
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Table 1 

The parameters for the GSRS force used in this work [33]. The second, the third, and 
the fourth columns are the interaction strengths v in MeV and the corresponding 
ranges 77 in fm for the tensor and tensor and LS channels, respectively. 
(See Eqs. (32) and (34).) 





Tensor 


LS 




3e 


30 


30 


vi (MeV) 


-7.5 


2.5 


-800 


m (fm) 


2.5 


2.5 


0.6 


V2 (MeV) 


-67.5 


20. 


800 


m (fm) 


1.2 


1.2 


0.4 


V3 (MeV) 


67.5 


-20. 




m (fm) 


0.447 


0.6 





Hartrcc-Fock (CPPHF) scheme, in which we perform both the parity and 
charge projections. We show the calculated results in Table 2, where we take 
xt — 1.5 and xte = 0.81, with which the binding energy of the alpha particle 
is reproduced in the CPPHF case. For comparison, we show the results for the 
HF calculation with the original Volkov No. 1 force {xt — 1.0 and xte — 1-0). 

In the simple HF case the energy from the tensor force (^t)*^^'^-* is zero. It 
means that the result of the HF calculation simply becomes a (Os)^ configura- 
tion and there is no p-state component. In this case the expectation value of the 
tensor force is zero identically, because the tensor force does not act between 
s-statcs. If we perform the parity-projection (PPHF), the energy contribu- 
tion from the tensor force becomes finite. The kinetic energy becomes larger 
because some component of the s-state is shifted up to p-states to gain the 
correlation caused by the tensor force. We perform then the charge projection 
(CPPHF) further. We see the contribution from the tensor force becomes much 
larger. It is reasonable because in the PPHF case only the rj^r^ component of 
Ti • T2 in the tensor force is active, while in the CPPHF case all the T1T2, 
T°T2 , and T1T2 parts of the tensor force are active. In fact, the energy from 
the tensor force in the CPPHF case is about three times large than that in the 
PPHF case. The results indicate that both the parity and charge projections 
are very important to treat the tensor force in a mean-field-type model. 

We note that the variation-after-projection scheme that we take here, is needed 
because even if we assume the mixing of parity and charge in the simple 

Hartree-Fock calculation, we cannot obtain the result with mixed symmetries. 
The parity mixed state contains odd parity states, which require large en- 
ergy jump across a major shell. In fact, the energy gap between Os-shell and 
Op-shell in the alpha particle is more than 20 MeV. Therefore, we need to per- 
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Table 2 

The results for the ground (0+) state of the alpha particle for various cases. HF de- 
notes the simple Hartree-Fock scheme. PPHF denotes the parity-projected Hartree- 
Fock scheme in which only the parity-projection is performed. CPPHF denotes 

the charge- and parity-projected Hartree-Fock scheme in which both the charge 
and parity projections are performed. The potential energy (('0)^^'^^), kinetic en- 
ergy ((T')(+'^)), total energy root-mean-square matter radius (Rm), and 
the probability of the p-state component (P(-)) are given in the table. {vc)^'^''^\ 
{vt)^^'''^\ {vls)^^'''^\ and (t'Coul)*-"'''^'* are the expectation values for the central, ten- 
sor, LS, and Coulomb potentials, respectively. The factors, xt and xte, denote the 
factors multiplied to the ti • T2-type tensor force and the central force. 





HF 


HF 


PPHF 


CPPHF 


XT 


1.0 


1.5 


1.5 


1.5 


XTE 


1.0 


0.81 


0.81 


0.81 


(i)c>(+'2) (MeV) 


-76.67 


-56.85 


-61.31 


-64.75 


{vt)^+'^^ (MeV) 


0.00 


0.00 


-10.91 


-30.59 


{VLS)^+'^^ (MeV) 


0.00 


0.00 


0.67 


1.91 


(^)coui)^+'') (MeV) 


0.83 


0.76 


0.78 


0.85 


(t))(+'2) (MeV) 


-75.84 


-56.10 


-70.76 


-92.58 


(f)(+'2) (McV) 


48.54 


39.98 


49.67 


64.39 


(MeV) 


-27.30 


-16.12 


-21.09 


-28.19 


Rm (fm) 


1.48 


1.63 


1.51 


1.39 


P(-) (%) 


0.0 


0.0 


7.6 


16.1 



form the parity projection before variation to remove the spurious component 
mixing into the intrinsic wave function ^^^'^ in (28) properly. In contrast to 
the parity mixing, the correlations by deformation and pairing can be treated 
quite nicely with the projection-after-variation scheme including a configura- 
tion mixing along the quadrupole degree of freedom [26,27,28,29], where the 
rotation symmetry and the particle-number symmetry are broken in a intrin- 
sic state. The correlations by deformation and pairing are in principle those 
within the same major shell and causes small energy loss and therefore the 
mixing of the angular momentum and the particle number occur even in the 
projection-after- variation scheme. 

We check now the dependence of the results on the strength of the tensor 
force by changing xt from 1.0 to 2.0. We determine xte so as to reproduce 
the binding energy of the alpha particle for each Xt- We show the results of the 
CPPHF scheme in Table 3. From the table we can see that the contribution of 
the energy from the tensor force becomes larger with xt- The kinetic energy 
becomes larger also. The probability of the p-state (P(-)) changes from 10.5% 
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Table 3 

The results for the ground (0"*") state of the alpha particle with the CPPHF scheme. 
We change xt from 1.0 to 2.0 and determine xte accordingly so as to reproduce the 
binding energy of the alpha particle. The potential energy 

kinetic energy 

((T)'^'^-*), total energy (E^^''^^), root-mean-square matter radius (i?m)) root-mean- 
square charge radius (Re), the probability of the p-state component (P(— )), and 
the probability of the 5 = 2 component {P{S = 2)) are shown. {vc)^'^'''^\ {vt)^'^'''^\ 
(?)ls)^"'"'^\ and {vcoui)^'^'^^ are the expectation values for the central, tensor, LS, 
and Coulomb potentials, respectively. 



XT 


1 


1.25 


1.5 


1.75 


2 


XTE 


0.93 


0.88 


0.81 


0.73 


0.64 


(uc)(+>2) (MeV) 


-73.60 


-70.13 


-64.75 


-58.34 


-50.69 


{v^)(+'^) (MeV) 


-12.26 


-20.23 


-30.59 


-43.86 


-60.41 


(t)Ls)(+'') (MeV) 


0.75 


1.26 


1.91 


2.72 


3.72 


(^Coui)(+'') (MeV) 


0.85 


0.85 


0.85 


0.86 


0.87 


(MeV) 


-84.27 


-88.26 


-92.58 


-98.62 


-106.51 


(f)(+;2) (MeV) 


55.81 


59.71 


64.39 


70.52 


78.19 


E(+;2) (MeV) 


-28.46 


-28.55 


-28.19 


-28.10 


-28.32 


Rm (fm) 


1.43 


1.41 


1.39 


1.36 


1.33 


Rc (fm) 


1.64 


1.62 


1.60 


1.58 


1.56 


P{-) (%) 


10.5 


13.3 


16.1 


18.8 


21.4 


P{S = 2) (%) 


2.9 


4.9 


7.3 


10.1 


13.3 



to 21.4%. In the last row in Table 3, we show the probability of the total spin 
S = 2 component P{S = 2), which is defined as, 

^(5=2)=„i:^ w^w^) ■ <^«' 

where 

^s=2iMs) = [[xi X X2]^'^ X [X3 X x.MSs- (37) 

P{S = 2) corresponds to the D-statc probability, because in the present calcu- 
lation the total wave function has the 0"^ spin-parity. P{S — 2) changes from 
2.9% to 13.3% when xt changes from 1.0 to 2.0, and is almost proportional to 
the tensor correlation energy (■Ot)*-^'^-'. Wc note that the D-state probability of 
the alpha particle in the recent exact-type calculations [35] is about 15%. The 
matter root-mean-square radius i?m decreases with xt- The root-mean-square 
charge radius Rc is calculated from the proton root-mean-square radius Rp 
as Rc — ^Rp + 0.64. This approximation for Rc corresponds to assuming the 
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charge radius of proton as 0.80 fm. The charge radii obtained here are shghtly 
smaller than the experimental value Rc = 1.676(8) fm [36]. The charge form 
factor has a dip around the momentum transfer squared = 15 fm~^ for 
the case without the tensor force, and with the inclusion of the tensor force 
the position of the dip moves towards smaller momentum transfer. The posi- 
tion of the dip is around = 11 fm~^ for ,xt = 1.5 in accordance with the 
experiment, but the amount of the second bump is somewhat underestimated. 



l.On 




Fig. 1. Single particle wave functions squared for the case with xt = 1-5 and 
xte = 0.81 as a function of the radial coordinate r (fm). The solid curve de- 
notes the positive-parity neutron component |</>;^o-iP = 0), the dotted curve the 
negative-parity neutron component {I = 1), the short-dashed curve the 

positive-parity proton component |</'ioiP (^ = 0); a-^d the long-dashed curve the 
negative-parity proton component l^^^ip (I = 1). For comparison, we show the 
harmonic oscillator (HO) wave functions with the oscillator length 1.37 fm for the 
Os state (dashed and dotted) and the Op state (dashed and double dotted). 

In Fig. 1, we show one of the intrinsic single-particle wave function (squared 
wave function) for the alpha particle with xt — 1.5 in the CPPHF method. 
The dominant component of this wave function is the s-state (positive-parity 
state) and neutron. Due to the charge and parity mixings, there appear the 
components of the p-state and neutron, s-state and proton, and p-state and 
proton. In the wave function the probability of the negative parity component 
is 16% and that of the proton component is 17%. There is one more intrinsic 
single-particle state (n = 2), which has the s-state-and-proton component as 
the dominant one. In the figure the harmonic oscillator (HO) wave functions 
with the oscillator length b = 1.37 fm for the Os and Op state are also plotted 
for comparison. The Os-state HO wave function almost corresponds to the 
calculated wave function in the simple HF scheme with the original Volkov 
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interaction (xt = 1-0, Xte = 1-0). The s-state and neutron component in 
the CPPHF method has almost the same width as the HO wave function but 
the magnitude is somehow reduced. This reduction is caused by the parity 
and charge mixings. The mixing of parity reduces the wave function around 
the origin. It is interesting to see that the negative parity components have 
narrower widths as compared to the Op-state HO wave function [37]. If we cal- 
culate the overlap between the p-state wave function in the CPPHF method 
with the Op-state HO wave function changing the oscillator length bop, the 
maximum overlap is achieved at feop — 0.81 fm. It means that the negative 
parity component in the CPPHF method is not a simple Op state. The mix- 
ing of a higher momentum component is necessary to make the p-state have 
such a compact distribution. This fact suggests that the tensor force induces 
higher momentum components, which are not included in a simple shell model 
configuration. The mixing of the higher-momentum component results in the 
increase of the kinetic energy as seen in Table 3. The correlation induced by 
the tensor force produces more attractive energy and therefore the mixing of 
the higher-momentum component is favorable in total. 

Table 4 

The results for the first 0~ state in the alpha particle. The potential energy {{v)^~'''^^), 
kinetic energy ((r)(~'^)), total energy {E^~'^^), root-mean-square matter radius 
(Rm), and the probability of the p-statc component (P(— )) are shown. {vc)^~''^\ 
{vt)^^'''^\ {vhs)^^'''^\ and («coui)*-^'^^ are the expectation values for the central, ten- 
sor, LS, and Coulomb potentials, respectively. 



XT 


1 


1.25 


1.5 


1.75 


2 


XTE 


0.93 


0.88 


0.81 


0.73 


0.64 


(*c)(-'2) (MeV) 


-29.05 


-29.73 


-29.36 


-28.16 


-26.03 




-4.57 


-8.68 


-14.00 


-20.49 


-28.18 




0.38 


0.73 


1.11 


1.52 


1.96 


(^)Coul)(-'') 


0.39 


0.48 


0.54 


0.58 


0.61 


(^)(-;2) 


-32.85 


-37.21 


-41.72 


-46.55 


-51.63 


(f)(-;2) 


30.15 


34.43 


38.99 


43.80 


48.80 


£;(-;2) 


-2.70 


-2.78 


-2.73 


-2.76 


-2.83 


Rm (fm) 


3.23 


2.74 


2.40 


2.19 


2.06 


P(-) (%) 


10.7 


15.1 


19.3 


22.8 


25.8 



It is an interesting subject to solve the 0~ state in the CPPHF method, since 
the 0~ state is the daughter state of the 0+ ground state in the parity- and 

charge-mixed intrinsic state. We calculate the 0~ state of the alpha particle 
with the same parameters as the 0^ state. The results are tabulated in Table 4. 
The attraction due to the tensor force increases with the tensor parameter xt, 
but the amount is less than a half of the case of the positive parity state. 
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The net binding energy comes out to be about 2.8 MeV. It is interesting to 
note, however, that the admixture of the opposite parity components P(— ) 
in the single-particle states is larger than the case of the O"*" state as seen 
in Table 3. This is related with the fact that -P(— ) ought to be finite for 
the 0~ state even for the case without the tensor force, while for the 0"*" 
case -P(— ) becomes zero in such From Table 3 and 4 the excitation 

energy from the ground state of the 0^ state in the CPPHF method is 
about 25.6 MeV in our calculation. The experimental value for the first 0~ 
state, which corresponds to the 0~ state calculated here, is Ej; — 21.010 MeV 
with the total width F = 0.840 MeV [38]. In the case of the alpha particle, 
the energy difference between the daughter states, the positive and negative 
parity states, is reasonably large, although the intrinsic single-particle states 
have large breaking of the parity and charge symmetries. 




r (fm) 

Fig. 2. Densities for the alpha particle with xx = 1-5 and xte = 0.81 as a function 
of the radial coordinate r (fm). The density of the 0^ state (the solid curve) and 
that of the 0~ state (the dotted curve) together with that of the 0"^ state with the 
original Volkov No. 1 force in the simple Hartree-Fock scheme (the dashed curve) 
are shown. 



Finally, we show the density distributions of the 0+ and 0^ states in Fig. 2. 
The density of the 0+ in the CPPHF scheme decreases around the center and 
increase in the middle region (around r = 0.8 fm) from that of the HF scheme. 
This change is induced by the admixture of the p-state component. In fact, 
the wave function squared for the p-state in Fig. 1 have the maximum around 
r = 0.8 fm. Due to the smaller binding, the density of the 0" state is much 
more spread over the space than the case of the 0+ state. 
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4 Summciry 



We have developed a mean-field framework with the projection method in 
order to treat the tensor force. To do this we mix parities and charges into a 
single-particle state, instead of using a single-particle state with a good parity 
and a definite charge number. We perform the projections of parity and charge 
from the intrinsic wave function, consisting of the parity- and charge-mixed 
single-particle states. We take the variation of the energy expectation value of 
the Hamiltonian with the projected total wave function, and obtain a Hartree- 
Fock-hke equation, the charge- and parity-projected Hartree-Fock (CPPHF) 
equation. 

We have applied the CPPHF equation to the alpha particle and shown that 
the CPPHF method is able to treat the strong tensor correlation. The varia- 
tion after projection is very important to obtain the tensor correlation in the 
present framework, because the mixing of parity and charge does not occur if 
the parity and charge projections are not performed before variation. 

We have obtained a large tensor correlation energy in the CPPHF method. The 
tensor correlation energy is about 31 MeV for xt = 1-5. To obtain such large 
tensor correlation energy, the charge mixing and projection on top of the parity 
mixing and projection is very important due to the isovector character of the 
one-pion-exchange interaction. We have found that if we perform the mixing 
of both charge and parity and their projections, the tensor correlation energy 
becomes about three times larger than in the case where only the parity mixing 
and projection is performed. In the result of the alpha particle, we see that the 
p-state component has the narrower width than the s-state component, and 
does not seem like a simple Op-state. It means that to obtain the correlation 
of the tensor force we need to mix higher-momentum components. This fact 
indicates the increase of the kinetic energy by the inclusion of the tensor 
correlation. In spite of the increase of the kinetic energy, the correlation energy 
from the tensor force produces more attractive energy and a net energy gain 
is achieved in total. The probability of the p-state component in the single- 
particle states becomes 10.5% to 21.4% with the parameters adopted here. 
The probability of the total spin S — 2 component, which corresponds to 
the D-state probability, changes 2.9% to 13.3% with the same parameters. As 
for the negative parity state (0^), the total binding energy becomes around 
2.8 MeV. We have found that the density for the 0+ state decreases around 
the center and increases in the middle-range region due to the mixing of the 
p-state component. The density distribution for the 0~ state is much broader 
than the ground state (0+). 

In this paper, we have pursued a mean-field (single-particle) framework with 
the projection method, in which the tensor force can be treated in the same 
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way as the central and LS forces. We have made such a framework by combin- 
ing the mixings of charge and parity in a single-particle state and the charge 
and parity projections on the total wave function. In the present study, we 
have apphed the CPPHF method only to the most simple closed-shell nucleus, 
the alpha particle. The apphcation of the CPPHF scheme to heavier closed- 
shell nuclei, for example, ^^C and ^^O, will show the role of the tensor force 
in the formation of the shell structure of nuclei. The extension of the present 
framework to the deformed case is also interesting and will give us the insight 
of the role of the tensor force on the nuclear deformation and the nuclear 
clustering. An application of the present method on the relativistic mean field 
framework is under progress [39]. 
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